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Abstract 

Within Bohm's interpretation of quantum mechanics particles follow 
"classical" trajectories that are determined by the full solution of the 
time dependent Schrodinger equation. If this interpretation is consistent 
it must be possible to determine the probability distribution at time t, 
p(x,t), from the probability distribution at time t — 0, p(x,0), by using 
these trajectories. In this paper it is shown that this is the case indeed. 



1 Introduction 



In 1952 David Bohm W proposed a new interpretation of quantum mechanics 
connecting "classical" trajectories to particles with a space-time probability dis- 
tribution given by the full solution of the time dependent Schrodinger equation. 
A detailed presentation of this interpretation may be found for instance in Q , 
whereas its historical contingency is discussed in ||. 

The main points of his interpretation, which are of concern here, may be 
described as follows (for simplicity the space dimension is taken to be 1 and 
units H = m = 1 are used). Let tp{x,t) be the solution of the time dependent 
Schrodinger equation for a certain system, that is 

i^(x,t) = Hi/>(x,t) (1) 

where H denotes the Hamilton operator of the system. Then one decomposes 
ip{x,t) in its absolute value and the corresponding phase which results in 

ip(x,t) = R(x,i)exp[iS(x,t)] (2) 

The phase S(x, t) determines the "classical" motion by requiring that the mo- 
mentum of a particle at position x at time t is given by 

dx dS(x,t) 

p= Tt=^x- (3) 

Solving (^) for a particle with initial position xq at time to attributes to this par- 
ticle a classical trajectory x(xq, t) subject to the initial condition xq — x(xq, to). 
Obviously one may use these trajectories to transport the probability distribu- 
tion 

p(x,t ) = |^(Mo)| 2 =R(x,t Q ) 2 (4) 

in time by moving the initial distribution at points xq to the points x(xo,t), 
which is suggested by Bohm's interpretation of quantum mechanics. Of course 
this transport in time cannot be the complete procedure since the normalization 
condition of the probability distribution must be maintained too. In the next 
section it will be shown how one can reconstruct the probability distribution 
p(x,t) from the initial one, p(x,tg), by using classical trajectories, the solutions 
of (||), and that the result agrees with the usual quantum mechanical one, namely 

P (x,t) = |V0M)| 2 (5) 
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2 Reconstruction of the Probability Distribu- 
tion 



Let the probability distribution at time t = be (t is chosen to be for 
simplicity) 

(6) 



p(x,0) = / d(x - x )f(x )dx = f(x) 



This is the continuous analogue of particles distributed at points xq with some 
probability f(xo)- According to Bohm's interpretation the points Xo move to 
the points x(xo,t) in time and therefore following this idea the probability dis- 
tribution at time t must be given by 



p(x,t) = / 8((x - x(x ,t))f(x )dx 



(7) 



where x(xq, t) is the solution of (g) with the appropriate initial condition. From 
(0) it is easy to obtain the wanted result because the integration can be done 
without difficulty by means of the 5- function. The zero of the d- function is given 
by xq(x, t), the inverse function to x(xo, t). This gives, using differentiation rules 
for inverse functions, 



p(x,t) 



dxo(x, t) 



and the final result is 



p(x,t) 



dx 



8xq(x, t) 



f(x (x,t)) 



0:i 



p(x (x,t),0) 



(8) 



(9) 



That this result, obtained by using Bohm's interpretation of quantum mechan- 
ics, agrees with the usual quantum mechanical result (||) can be shown as follows. 
Together with the current density 



dx(xo, t) 
dt 



S(x - x(x ,t))f(x ,t)dx 



p(x,t) as defined in (^) fulfills the continuity equation 

dp{x,t) dj(x,t) 



= 



(10) 



(11) 



dt dx 

But j(x, t), due to the momentum definition can also be written in the form 



i(M) = 



dS(x,t) dS{x,t) 

6[x - x{xo,t))f{x ,t)dx = — p(x,t) 



dx 



dx 



(12) 
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Thus p(x, t) fulfills the same continuity equation as p(x, t) defined in 
time t = both agree, if f(x) is chosen to be 



At 



/(x)H^,0)| 2 (13) 

The uniqueness of the solution of the continuity equation implies then that 

p(x,t) = p(x,t) if p(x,0) = p(x,0) (14) 

This shows that the reconstruction of the probability distribution at time t — 
following Bohm's ideas agrees with the usual quantum mechanical result and 
the explicit expression is given by <^). But the imagination that the probability 
distribution at time t ist just obtained by the flow of particles from the initial 
distribution following the trajectories x(xo,t) is in general wrong due to the 



facte 



dxo(x,t) 
dx 



in (M). This factor could also be interpreted as a normalization 



factor since due to it the normalization of the wave function is valid for all times 

dxo (x, t) 



1= / |V(^ o ,0)r rfx = / |^(sso(;M),0)| 



dx 



dx 



(15) 



The next section will show the role of this factor within three different examples 
explicitly. 



3 Examples 

3.1 The Coherent State of the Harmonic Oscillator 

The normalized coherent state with amplitude d is given by (w = 1) 



1 

— exp 

7T4 



t 



sin 2i + dx sin t 



— \x — d cos t] 2 
2 L J 



2 4 

Which gives for the absolute value R(x,t) and the phase factor S(x,t) 



R(x,t) 



S(x,t) 



1 



■ exp 



t 



— (x — dcost) 2 



■ sin 2t — dx sin t 



2 4 

For the definition of the momentum one obtains therefore 

p = — — = —a sin t 

ox 

resulting in the equation of motion 

dx 



dt 



= —dsmt 



(16) 

(17) 
(18) 

(19) 
(20) 
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with the solution 

x(xo,t) = d (cost — 1) + xq 
From this the inverse function is given by 

Xq(x, t) = x — d (cost — 1) 

and therefore 



dx (x, t) 



dx 



1 



(21) 
(22) 
(23) 



as expected since a coherent state does not change its shape in time. Obviously 
the reconstructed probability distribution agrees with the quantum mechanical 
result as it has to be 



1 r , 

p(x,0) = R(x,0) = — cxp — (x — dy 

7T2 L 

p(x (x,t),0) = — exp <^ - [xo(x,t) - df \ — R(x,ty 

7T2 I- - 1 



(24) 



3.2 The Wave Packet for a Moving Free Particle 

The normalized wave function for a free particle with a Gaussian distribution 
around x — at time t = and for which the peak of the probability distribution 
moves with constant velocity (here taken to be one) in time is given by 



ip(x,t) 



1 



7T 4 y/ 1 + it 



exp 



x 2 - 2ix + it 
2(1 + it) 



(25) 



From this the absolute value and the phase of the wave function are given by 



R(x,t) 



i/(t 2 + l)TT 



exp 



1 (t - x y 

2 t 2 + l 



„. . 1 1 - it tx 2 + 2x - t 

with the corresponding "classical" momentum 

_ dS _ tx + 1 

p ~ ~dx~ ~ TTt2 



The equation of motion 



has the solution 



dx tx + 1 
~dt = 1 + t 2 

x = t + x \Jl + t 2 



(26) 
(27) 

(28) 

(29) 
(30) 
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whose inverse is given by 



and 



x = 



x — t 



dx n (x, t) 



da 



VT+7< 



(31) 



(32) 



which is just the right factor to correct the normalization for the broadening of 
the wave packet in time. Again one obtains complete agreement with Bohm's 
interpretation 



p(x,0) 



R(x,0) 2 = cxp(-a; 2 ) 

V 7 *" 



8xq(x, t) 



exp 



(x - if 



l + t 2 



dx 
R(x,tf 



p(x o (x,t),0) (33) 



One may be tempted to interpret the factor g*' to be responsible only for 
the normalization of the wave function without any relation to its shape. That 
this imagination is completely wrong will be shown in the last example. 



3.3 The Harmonic Oscillator Again 

Instead of a coherent state a superposition of the ground state with the first 
excited state of the harmonic oscillator is considered. The time dependent 
solution for this case is given by 



ip(x,t) 



2 1 



3^/4 



exp 



x 2 + it 



[1 + a; exp (—it)] 



and therefore 



/2 1 / x 2 \ i 

R(x, t) = \ — —rjj exp — — v 1 + 2x COS t + ; 
V 3 it i \ 2 / 

t 1 1 + x exp (-it) 
b(x,t) = — — + — in — ; — — 

v ' ' 2 2i l + xcxp(it) 



(34) 



(35) 



(36) 



Proceeding in the usual way one has therefore to solve the differential equation 

dx sin t 

dt 1 + 2x cos t + x 2 

which has the implicit solution 



(37) 



^\/7rerf(a;) — exp (—x 2 ) cost — ^xcxp (—x 2 ) 

3 1 

-V^erf (x ) — exp (— x 2 ,) cost — -x exp (— x%) 



(38) 
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Figure 1: Comparison of the exact result R(x,ir) and the result obtained by 
using "classical" trajectories but without the "normalization" factor 



One can numerically extract from this implicit solution x(xo,t) or the inverse 
xq(x, t) and compare the exact result p(x, t) with the result obtained from using 



dxp(x,t) 
dx 



"classical" trajectories but without the factor 
only. This is shown in Fig. 1 for the time t = 
roots are plotted, that is R(x, t) and R(xo(x, t), 0)) 

is obvious, but full agreement is achieved if the factor 
account. 



, that is with p(xo(x, t), 0) 

(for simplicity the square 
The disagreement in shape 

is taken into 



dxo(x,t) 
dx 



4 Conclusion 

It has been shown that Bohm's interpretation of quantum mechanics gives the 
correct answer for the probability distribution at later times when applied to 
the probability distribution of a system at time t = t . Of course no new results 
can be obtained in this way, since first one has to know the wave function of 
the system for all times. But the main point is that the interpretation in itself 
does not lead to any contradictions if one takes it literally. Nevertheless one 



G 



has to be very cautious in applying this interpretation. As is shown by the last 
example of the superposition of two harmonic oscillator eigcnstatcs, one cannot 
view the particle flow as just transferring the initial proabability distribution 
to a final one, which has to be normalized only, but agrees in shape with the 

may change 



dxo(x.t) 



dx 



result obtained in this way. On the contrary, the factor 

this shape quite drastically in general and makes it very difficult to visualize 
what is going on. 
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